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ABSTRACT 

We have observed the radio nebula surrounding the Galactic LBV candidate 
G79.29+0.46 with the EVLA at 6 cm. These new radio observations allow a 
morphological comparison between the radio emission, which traces the ionized 
gas component, and the mid-IR emission, a tracer of the dust component. The 
IRAC (8/im) and MIPS {2Afim and 70 /im) images have been reprocessed and 
compared with the EVLA map. We confirm the presence of a second shell at 
24/im and also provide evidence for its detection at 70 /im. The differences be- 
tween the spatial morphology of the radio and mid-IR maps indicate the existence 
of two dust populations, the cooler one emitting mostly at longer wavelengths. 
Analysis of the two dusty, nested shells have provided us with an estimate of the 
characteristic timescales for shell ejection, providing important constraints for 
stellar evolutionary models. 

Finer details of the ionized gas distribution can be appreciated thanks to the 
improved quality of the new 6 cm image, most notably the highly structured 
texture of the nebula. Evidence of interaction between the nebula and the sur- 
rounding interstellar medium can be seen in the radio map, including brighter 
features that delineate regions where the shell structure is locally modified. In 
particular, the brighter filaments in the south-west region appear to frame the 
shocked southwestern clump reported from CO observations. 

Subject headings: circumstellar matter — infrared: stars — stars: early-type — stars: 
individual (G79.29+0.46) — stars: winds, outflows 
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Introduction 



Massive stars play a fundamental role in the evolution of galaxies. They are major 
contributors to the interstellar UV radiation and, via their strong stellar winds, provide 
enrichment of processed material (gas and dust) and mechanical energy to the interstellar 
medium. Despite their importance, the details of post-MS evolution of massive stars are 
still poorly understood. Recent evolutionary models suggest that Luminous Blue Variables 
(LBVs) and related transition objects may play a key role in the massive sta r evolution, 



repre senting a crucial phase during which a star loses most of its H envelope (JLamers et al 
2OOII ). More recently, it has also been pointed out t hat LBVs might be direct Supernovae 



progenitors (JKotak Sz Vink 



2006 



Smith et al. 



20081 ). enhancing their importance in the 



framework of stellar evolution 

LBVs are luminous (intrinsically bright, L ~ 10^ L©) objects, exhibiting different kinds 
of photometric and spectroscopic variabilities. They are massive (M ~ 20 — 120 M©), 
characterized by intense mass-loss rates (10~^ — 10~^MQyr~^), which can also occur in the 
form of eruptive events. Although eruptive events have been witnessed very rarely (e.g., 
r] Car and P Cyg) the presence of extended, dusty circumstellar neb ulae around LBVs 



(LBVNe) suggests that they are a common aspect of LBV behavior ( IWeisI l2008l ) . There 
are, however, many aspects of LBV phenomenon that are not completely understood. 
Among these are the total mass lost during the LBV phase (a key parameter necessary to 
test evolutionary models), the origin and shaping of the LBVNe, and how the mass- loss 
behavior (single versus multiple events, bursts) is related to the physical parameters of the 
central object. 

The mass-loss archeology of the central object can be recovered from an analysis of 
its associated nebula. A successful approach is based on a synergistic use of different 
techniques, at different wavelengths, that allows one to analyze the several emitting 
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components coexisting in the nebula. In particular, a detailed comparison of mid-IR and 
radio maps, with comparable spatial resolution, has provided estimates of both ionized gas 
and dust masses and allowed us to sort out morphological difference s in the maps whic h 



can be associated with mass-loss behaviour during the LBV phase ( iBuemi et al. 



Umana et al. 



2OIOI). 



2010 : 



1.1. The nebula surrounding G79. 29+0.46 



G79.29+0.46 is considered a LBV candidate because its observed properties to date 
do not meet the requirements of spectral and photometric variability to be accepted as a 
bona-fide member. However, H^ variation s, suggestive of ma ss-loss variability during S-Dor 



variations, have recently been reported by 



Vink et al 



(120081 ). 



T he highly symrn e tric r inglike structure surrounding G79.29+0.46 was first pointed 



out by 



Wendker et al. 



determined by 



(11991). T 



Higgs et al. 



le thermal nature of the continuum radio emission was 



(J1994J ). who suggested that the ringlike nebula is an ionized 



shell of swept-up in terstellar material. By examining the IRAS high-resolution images. 



Waters et al. 



(119961 ) concluded that the ringlike nebula is a detached shell formed during 
an epoch of high mass loss (~ 5 x 10~^ MQyr^^) followed by a quieter period. This scenario 
is consistent with G79.29+0.46 being a LBV, where different mass loss events may have 



occurred in the recent past. 



More recently, the Spitzer Space Telescope has provided high 
IRAC and MIPS maps (The Cygnus-X Spitzer legacy p rogram, 



provides a better understanding of the dust properties (IKraemer et a 



sensitivity. 



Hora et al. 



the detection of CO millimetre emission in G79.29+0.46 ( iRizzo et al 



ligh- resolution 



2OIOI) which 



20101) . Moreover, 



20081 ) demonstrates 



that another component, consisting of molecular gas, is present in the surroundings of this 
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star and should be included in the budget of t 



In their CO maps of G79.29+0.46, 



l e tota l mass lost from the central object. 



RizzoetaL 



(|2008[ ) identified components of warm and 



dense molecular gas whose morphology closely resembles that of the extended IRAS nebula. 
They also reported the presence of a shock front, interpreted as a natural consequence of 
different wind regimes during the central object's evolution. 

Despite the wealth of new mid-IR and mm observations of the sour ces, our knowledge 



Higgs et al. 



(119941 ) 



of the radio emission has been limited to the 1988 VLA data reported by 

and on the 1400 and 350 MHz Westerbork observations, carried out between 1996 and 

1997 and reported by 



Setia Gunawan et al 



(120031 ). In this paper, we present new EVLA 



observations with sufficient dynamical range, sensitivity, and angular resolution to provide 
a good match to the Spitzer images. The EVLA observations allow, for the first time, a 
detailed morphological comparison with other maps tracing the different components of the 
nebula. 

2. Observations and Data Reduction 
2.1. EVLA observations 



Radio observations of G79.29+0.46 and its associated nebula were carried out with the 
EVLAI^ on 2010 June 1 and 11 and 2010 December 1 and 5. The source was observed at two 
frequencies (1.4 and 5 GHz), with a total bandpass of 256 MHz and in two different array 
configurations (D and C). For each configuration and frequency, the same observing scheme 
was followed; namely, a 5 minute on-source scan preceded and followed by a 1 minute scan 
on the phase calibrator, for a total of 75 minutes on-source integration time. The source 



^The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 
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J0137+3309 was used to fix the absolute amplitude scale and to correct for the bandpass. 
Table 1 contains a summary of the observing details. The data reduction was carried out 
within the Common Astronomical Software Applications (CASA) package, version 3.0.2. 
At each frequency, the data from each configuration were first independently calibrated and 
then combined into a single UV data set. The imaging process was performed by setting 
the Briggs robust parameter equal to 0, a compromise between uniform weighting of the 
baseline for highest angular resolution and natural weighting for highest sensitivity. We also 
used a multi-scale CLE ANing algorithm, in tended for high resolution image but sensitive 



to extended structures f lBrogan et al.ll2006l ). resulting in a single image with a rms noise of 
0.07 mJy beam~^ and a synthetic beam of 4'.'56x3'.'09 for the 5 GHz observations, and a rms 
noise of 0.5 mJybeam"^ and a synthetic beam of 18'.'5x 18'.'6 for the 1.4 GHz observations. 



2.2. Spitzer data 



2.2.1. IR AC data 



Infrared imaging of the field including G79.29- 



8.0 /im wit 



1 the InfraRed Arr ay Camera (IRAC) (JFazio et al.ll2004j ) on the Spitzer Space 



0.46 was perform ed at 3.6, 4.5, 5.8, and 



Telescope ( iWerner et al 



2004J ). All available data for G79.29+0.46 from the cryogenic 
Spitzer mission archive were used, including AORIDs 6050560, 17330688, 27106560, 
and 27107584. The observations used the 12 and 30 second HDR modes, which obtain 
integrations with frame times of 0.6 & 12 seconds and 1.2 & 30 seconds, respectively. 
The Basic Calibrated Data (BCD) were retrieved from the archive (pipeline version 
S18.18), and additional interactive processing was performed using the imcleaiu tools 



^http://irsa.ipac. caltech.edu/data/SPITZER/docs/dataanalysistools/tools/ 
contributed/irac/imclean/ 
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to remove image artifacts from bri ght sources, including muxbleed, banding, and column 
pulldown e ffects (U 



Hora et al. 



IRACproc ( ISchuster et al. 



20041 ) ■ Mosaics at each wavelength were constructed using 



20061 ). combining the data which was taken at different spacecraft 



orientations and resampling to a final image pixel size of 0'.'6. Outliers due to cosmic rays 
or instrument artifacts not previously flagged are removed in this process. The images were 
combined using a weighting based on the integration time of the frames. In pixels that are 
saturated in the longer exposures, only the shorter frames are used in the final image. 



2.2.2. MIPS data 

The MIPS data was taken as part of the Cygnus-X Spitzer Legacy program (PID 
40184, PI J. Hora) in the fast scanning mode with a cross scan step of 148" to fill the 70 ^m 
array, resulting in an integration time of 15.7 seconds per pixel on a sin gle scan leg at bot h 



24 and 70 iJ,m. The 24 u m data were reproce ssed using the MIPSGAL ( ICarey et al. 



data pro cessing pipeline 



pipeline ( jPaladini et al. 



([Mizuno et al. 



20091) 



20081 ). The 70 /xm data also used the MIPSGAL 



201ll ) with the exception that the non-linearity correcti on at 70 um 



was done per pixel using the prescription defined by the SINGS Legacy team ( iDale et al 
20071 ). based on the behavior of the MIPS 70 fim calibrators. 



3. The radio nebula associated with G79.29+0.46 

As we intend to compare the spatial distribution of the ionized gas component, traced 
by the radio observations, with the morphology of the dust component, traced by the mid-IR 
observations, in the following we will analyze and discuss only the 5 GHz dataset, as that 
image provides details with a spatial resolution comparable to that of Spitzer observations. 
The morphology of the radio nebula is evident in our 6 cm EVLA multi-configuration 



image, whic h reveals a well def ined shell-like structure whose overall shape was previously 
reported by 



Higgs et al. 



( Il994l ). However, the improved quality of our image, shown in 



Figure [H allows us to discern finer details of the ionized gas distribution, most notably 
the highly structured texture of the nebula. There is evidence for interaction between 
G79.29-I-0.46 and the surrounding interstellar medium: i.e., the bright frontal structures 
in the northeast and the southwest regions of the nebula that delineate regions where the 
shell structure is locally modified. In particular, the brighter filaments i n the south-west 



region appears to frame the shocked southwestern clump observed in CO ( JRizzo et al. 



2008 



Jimenez-Esteban et al. 



2010h . 



The central object is well-detected, as is the (probably extragalactic) background 
object to the northwest. The position, flux density and angular sizes of these components 
have been derived by fitting two dimensional Gaussian brightness distributions to the 
map. We obtained a flux density of 1.51 ± 0.08 and 1.88 ± 0.09 mJy for the central and 
background components, respectively. The uncertainty associated with the flux density 
estimation is given by: 

a = ^{RMS,,,Y + (a,,i^,)2 (1) 



where RMStot is the rms noise in the map, and Ucai is the systematic error due to the flux 
calibrator (typically on the order of 3%). 



by 



Within the error s, the derived flux densities are in agreement with those determined 



Higgs et al. 



(J1994J ). who derived a spectral index of 1.39 ± 0.14 between 5 and 8.4 GHz. 
They pointed out that, even though this spectral index is steeper than the canonical 0.6, it 
is still consistent with mass-loss from a stellar object. 



Assuming that the central radio source is related to the ste 



can derive its current-day mass-loss from the standard formula (jPanagia fc FeUi 



lar wind from the LBV, we 



19751): 



M = 6.7 X lQ-\^Fl'^D%l{u X (7ff)-o-5 M^yr-i 



(2) 
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where full ionization and cosmic abundances have been assumed, Fj, is the observed radio 
flux density, in mJy, and Voo is the terminal velocity of the win d in kms~^. The fr ee-free 



y3/2 



Gaunt factor gg is approximated by gg = 9.77(1 + 0.13 log ^^ — ) ( JLeitherer fc Robertlll99ll ). 



From the radio fl ux density o 
V ^ 110 km s~ (Voors et al. 



Jimenez-Esteban et al 



Vink et al. 



jserv ed at 5 GHz, assuming as stellar wind velocity a value of 
2000l ). a wind temperature of lO'^K and a distance of 1.7 kpc 



20101), we derive a mass loss rate of M = 5 x 10^"^ Mqjt^^ . 
(l2008l ) pointed out the possibility that G79.29+0.46 is associated with the 
nearby DR15 region rather than with Cyg 0B2. If this is the case, G79.29-I-0.46 is likely to 
be closer, at about Ikpc. Since mass- loss scales with the distance as D'^/^, this would imply 
a reduced mass-loss of a factor 0.45. These values are smaller, but still consi stent, with 



previ ous evaluations obtained by other authors using different techniques (e.g., 



Waters et al. 



19961 ) , but it is at least an order of magnitude sm aller than the curren t-day mass- loss rates 



derived from radio measurements for other LBVs (lUmana et al.il2010l ) 



Ionized component versus dust components 



The 8;um IRAC and 24 and 70 fim MIPS images of G79. 2 9-1-0.46 are shown in 



Fig ure Ek- Previous versions of t hese images were presented by 



Kraemer et al. 



and 



Jimenez-Esteban et al 



f l2010h 



( 120101 ) ■ who both pointed out the presence of a second larger 
radius shell in the 24 fim image. Our reprocessed maps confirm the existence of the second 
24 //m shell and, thanks to our improved reduction, also provide a hint of its presence 
at 70 /xm. Both the 24 and 70 /xm maps show the same overall distribution of the dust, 
consistent with at least two nested dusty shells surrounding the LBV. On the contrary, in 
the 8 /xm map where a major contribution from warm dust is expected, the dust appears 
to be more concentrated in two south-west and south-east regions. The new EVLA radio 
observations allow, for the first time, a morphological comparison between the ionized gas 
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and the dust. In Figure [2)3, the same maps are shown but now with the 6 cm EVLA map 
superimposed using white contours. Once again there is a difference between the IRAC and 
MIPS maps: while the nebular emission at both 24 /im and 70 /xm is more extended than 
the ionized gas (radio nebula), the 8/im emission appears well-contained within the ionized 
part of the nebula. Moreover, while the radio emission from the nebula is incomplete in 
the north region, the emission at 24 /im and 70 /im has a more uniform distribution. The 
apparent lack of dust in the southern regions in both 8 /xm and 24 /xm are perhaps due to 
an absorption effect caused by the infrared dark cloud (IRDC) which may be located in 
front of the the southern region of the nebula. 

To better visualize the spatial distribution of ionized gas relative to the dust, surface 
brightness profiles have been extracted along cuts through the nebula from the 6 cm, 24 /xm 
and 70 /im maps. However, a direct comparison can be performed only between the 6 cm 
and the 24 /im maps, as they share a comparable angular resolution, ~ 5" for the 6 cm 
EVLA versus ~ 6" for the 24 /im MIPS map. A superposition of the 6 cm EVLA map 
on the 24 /im map is shown in Figure [3l We extracted eighteen cuts from each map and 
determined an azimuthally-averaged source profile, shown in (Figure H]). The radio shell 
has a smaller thickness and a sharp decreasing trend in the outer part of the nebula. The 
dust nebula is more extended and shows a smoother distribution in the outer regions. The 
second shell is clearly evident in the 24 /im profile, at a distance of about ±200" from the 
center. At the same distance, less defined peaks are visible in the 70 /im profile, smoothed 
out by the coarser spatial resolution. 

5. Discussion and Conclusions 

In this letter we present new radio EVLA observations of the nebula surrounding the 
LBV G79.29+0.46. In particular, we have analyzed the 6 cm map and, for the first time. 
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compared the radio free-free, which traces the spatial distribution of the ionized gas, with 
maps recently obtained in the mid-IR, which traces the spatial distribution of the dust 
component. 

The dust distribution observed at 8 /im appears different than the distribution at 24 
and 70 /zm. The warmer dust, traced by the 8/im emission, is more concentrated in the 
southwest and southeast regions. Moreover, when compared with the radio nebula, the 24 
and 70 /im emission appears more extended, while the 8 /im emission is well-contained inside 
the radio nebula. This is consistent with the existence of two dust components: a cooler 
component traced by the 24 and 70 /im emission, and a warmer dust component, probably 
consisting of smaller grains, traced by the 8/im emission. This warmer component, with 
morphological properties quite different from those obser ved at longer wavelengt h s, cou ld 



be related to the PAH emission that has been claimed by 



Jimenez-Esteban et al 



f!2010h 



An improved MIPS 70 /im map, constructed using the most updated pipeline from 
MIPSGAL team, detects the presence of a second shell very similar to that seen at 24 /im. 
Nebular lines and dust continuum are the main possible contributors to the 24 /xm emission 
from objects embedded in a dusty nebula with a hot central component. The possibility of 
such a combination, with one kind of emission being predominant, has been suggested by 



Mizuno et al 



( I2OIOI ) to explain the 24 /im emission detected in the MIPSGAL Bubbles. In 
the case of G79. 29-1-0. 46, the presence of similar shells at 24 /im and at 70 //m, together with 
the lack of any prominent emission line falling within the response c urve of the MIPS 24 urn 



band, as evident from the low-resolution IRS spectrum of the shell (jJimenez-Esteban et al. 



2OIOI ). leaves very little doubt that the shell emission is entirely due to thermal dust 



emission. 



The fact that there are at least two nested shells provides strong constraints on the 
the origin of the nebula, as it is difficult to explain the presence of multiple shells in the 
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hypot hesis of a nebula consisting of swept-up ISM material, as suggested by iHiggs et al. 



( Il994j ). It is evident that the dusty shells consist of material ejected by the central object 



in different mass-loss episodes. 

From the 24 /xm emission profile, we have determined that the inner shell peaks at 
100" and the second one at 200" from the central object. Such dust emission peaks can 



be re 



ated to the epoch when enh anced mass-loss took place. Assuming a distance of 1.7 



kpc f lJimenez-Esteban et al. 



20101), this corresponds to a linear distance of 0.82 pc and 



; .64 PC, respectively . This implies, assuming a shell expansion velocity of ~ 30 km sec 



( IWaters et al 



19961), that the two mass- loss episodes occurred 2.7 x 10 and 5.4 x 10 years 



ago, o r 1.6 X 10^ and 3.1 x 10^ years ago, if a distance of 1 kpc, as suggested by 



Vink et al. 



(120081 ). is used. Our results point out that mass- loss can occur in different episodes, for 
which we derived characteristics timescales, and pose quite strong constraints that stellar 
evolution models must take into account. 

The direct comparison of the 6 cm and 24 /im maps (that share comparable spatial 
resolution) indicates that only part of the inner, brighter nebula is ionized, (i.e. the nebula 
is ionization bounded). The radio nebula shows a sharper profile, depicting regions where 
probable interaction between the second and the first dusty shells is taking place. This 
is evident in several regions, most notably in the northeast and in the southwest part 
of the nebula, where the overall quite regular shell-like morphology is locally disturbed. 
The mid-IR dust emission depends linearly on density and as a modified blackbody on 
temperature; the mid-IR emission is therefore most sensitive to dust temperature. On the 
other hand, assuming an uniform temperature, the thermal radio emission depends on the 
square of the ionized gas density. This different dependence on density means that the 
radio emission will emphasize features that have the largest density. This is consistent 
with the spatial coincidence of the bright radio features with the CO emission reported 



by 



RizzoetaL 



(120081 ) and 
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Jimenez-Esteban et al. 



(J2010|). In particular, the high critical 



density of the J = 3 — 2 CO line [n^ ~ 10^ cm~^) indicates that the CO emission traces 
higher density regions, confirming that the brighter radio features delineate interactions 
between the shells where shocks can occur. 



This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. Support for this work was provided by NASA through an 
award issued by JPL/Caltech. 
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Table 1: Observational summary 



Observing EVLA Frequency Integration time Phase Cal 
dates Conf. GHz fmin) 



2010 June 1 


D 


1.4 


90 


J2052+3635 


2010 June 11 


D 


4.9 


90 


J2015+3710 


2010 Dec. 1 


C 


1.4 


90 


J2015+3710 


2010 Dec. 5 


C 


4.9 


90 


J2015+3710 



-15- 



A = 6 cm 



m 

c6 



150 



100 



50 







-50 



-100 



-150 




-150 -100 -50 



50 
asec 



N 

▲ 



E M- 



100 



150 



Fig. 1. — The 6 cm EVLA map of G79.29+0.46 obtained by combining the data from two 
different array configurations, C and D. The field of view is 2.5' x 2.5' centered on the LBV 
position. 
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Fig. 2.— Panel a: The 8/im IRAC (right), the 24 /xm (center) and the 70 /im (left) MIPS 
maps of G79.29+0.46. Panel b): the same as panel a, but the EVLA 6 cm map (white 
contours) superimposed. To better point out the major radio features, only contours starting 
from 10% of flux density of the central radio sources have been drawn. For each map, the 
field of view is 4.2' x 4.2' centered on the LBV position. 
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Fig. 3. — The 6 cm EVLA map of G79.29+0.46 (grey) superimposed to the MIPS 24 /im 
map (red). The field of view is 3.5' x 3.5' centered on the LBV position. 
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Fig. 4. — Averaged profiles through the nebula at 6 cm (dotted line), 24 /im (continuous 
line) and 70 /im (dotted slash line). Each profile has been obtained from 18 individual cuts 
in the corresponding map and has been shifted vertically by an arbitrary quantity for easier 
comparison. 
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